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ABSTRACT 

I suggest that the beaming factor in bright ULXs varies as b oc m~ 2 , where fa is the 
Eddington ratio for accretion. This is required by the observed universal L so ft oc T -4 
relation between soft-excess luminosity and temperature, and is reasonable on general 
physical grounds. The beam scaling means that all observable properties of bright 
ULXs depend essentially only on the Eddington ratio m, and that these systems vary 
mainly because the beaming is sensitive to the Eddington ratio. This suggests that 
bright ULXs are stellar-mass systems accreting at Eddington ratios of order 10 - 
30, with beaming factors 6>0.1. Lower-luminosity ULXs follow bolometric (not soft- 
excess) L ~ T 4 correlations and probably represent sitfe-Eddington accretion on to 
black holes with masses ~ 10M©. High-mass X-ray binaries containing black holes or 
neutron stars and undergoing rapid thermal- or nuclear-timescale mass transfer are 
excellent candidates for explaining both types. If the 6 am" 2 scaling for bright ULXs 
can be extrapolated to the Eddington ratios found in SS433, some objects currently 
identified as AGN at modest redshifts might actually be ULXs ('pseudoblazars'). This 
may explain cases where the active source does not coincide with the centre of the 
host galaxy. 
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1 INTRODUCTION 

There are currently two models proposed for ultralumi- 
nous X-ray sources (ULXs). In one they are identified as 
intermediate-mass black holes (IMBH) accreting at rates be- 
low their Eddington limits. In the alternative model, ULXs 
represent a very bright and unusual phase of X-ray bi- 
nary evolution, in which the compact object is fed mass 
at a rate M well above the usual Eddington value Me- In 
the picture proposed by Shakura & Sunyaev (1973) radi- 
ation pressure becomes important at the spherization ra- 
dius R sph ~ 27(M/Me)Rs/4:, where R s = 2GM 1 /c 2 is the 
Schwarzschild radius of the accreting black hole of mass M\ 
(Shakura & Sunyaev, 1973; see also Begelman et al., 2006; 
Poutanen et al, 2007). Within this radius the disc remains 
close to the local radiation pressure limit. Matter is there- 
fore blown away so that the accretion rate decreases with 
disc radius as M(R) ~ M(R/R sph ) ~ Me(R/R s )- As the 
disc wind has the local escape velocity at each radius, we 
see from mass conservation that the wind is dense near R sp h 
and tenuous near the inner disc edge, and there is a vacuum 
funnel along the central disc axis through which radiation 
escapes. 



In this model the large apparent X-ray luminosity 
Lx = 10 40 I/4o erg s _1 results from two effects of super- 
Eddington accretion (Begelman et al., 2006; Poutanen et 
al, 2000). First, the bolometric luminosity is larger than 
the usual Eddington limit by a factor ~ 1 + ln(M / Me), 
which can be of order 5 — 10 for the high mass transfer rates 
encountered at various stages of the evolution of a compact 
stellar-mass binary. Second, the luminosity of a ULX is col- 
limated by a beaming factor b via scattering off the walls of 
the central funnel. (Note that here and thoughout this pa- 
per, 'beaming' simply means geometrical collimation, and 
not relativistic beaming.) These conditions could occur in a 
state of high mass transfer (cf King, 2001, Rappaport et al., 
2005) or conceivably a transient outburst (King, 2002). In 
this picture one would expect on physical grounds that the 
Eddington ratio rh = M /Me should determine the beaming 
factor b. However current modelling has not yet derived this 
connection, allowing a spurious extra degree of freedom in 
comparing this picture with observations. 

A clue here comes from the fact that bright ULXs 
have spectra consisting of a power law plus a soft (kT ~ 
0.1 — 0.3 keV) excess which can be modelled as a black- 
body. This is usually taken as a multicolour disc with the 
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maximum disc temperature as the reference value, but the 
fitted temperature is not very different if the blackbody is 
assumed uniform. Feng & Kaaret (2007) show that the lumi- 
nosity L so f t and temperature T of the blackbody component 
vary as L so f t °c T~~ n with n = —3.1 ± 0.5 in the ultralu- 
minous source NGC 1313 X-2. Kajava & Poutanen (2008, 
hereafter KP)) have recently extended this result to a sam- 
ple of nine ULXs (including NGC 1313 X-2) which have a 
power law continuum with a soft excess. These are essen- 
tially all the sources with inferred luminosities permanently 
above ~ 3 x 10 39 erg s~\ Strikingly, KP find that all of 
these soft -excess objects cluster around the line 



L soCt = 7 x 10 40 T o ;ikcV erg 8 " 1 



(1) 



at all epochs (see the right-hand panel of their Figure 3). 
Here Tb.ikev is the temperature in units of 0.1 keV. KP 
caution that the agreement for the coolest and brightest may 
be affected by an incorrect subtraction of the hard emission 
component, but the overall trend ([1} is clear. 

KP also identify a distinct class of 'non-power-law' (or 
thermal) type ULXs whose medium energy spectra are fit- 
ted by harder multicolour disc blackbodies with reference 
temperatures feT me di um ~ 0.5 — 2 keV rather than power 
laws plus a soft excess. These systems all have inferred lu- 
minosities permanently below ~ 10 39 erg s _1 . They do not 
obey {TJ, but instead follow individual luminosity - tem- 



perature correlations L n 



, just like standard 



(non-ULX) black hole binaries (Gierlinski & Done 2004). 

At first sight, as KP remark, the correlation (JXJ for the 
soft excesses of bright, power-law ULXs seems counterintu- 
itive, as one might expect a blackbody to vary as L so ft tx T 4 . 
However this assumes that the characteristic radius R of the 
blackbody remains fixed as other parameters vary, and in- 
deed that the inferred L so ft is not affected by beaming, which 
could itself also vary systematically. 

I shall show here that in the picture of ULXs as super- 
Eddington accretors suggested by Begelman et al (2006) and 
Poutanen et al. (2007), the correlation L ao f t cx T~ 4 is actu- 
ally expected, and results from a tight relation between the 
beaming factor b and the Eddington ratio m of the form 
b ~ rh~ 2 . Using the observed relation JTJ we find b and 
mi = Mi/Mq as functions of the Eddington ratio m for a 
given inferred disc luminosity. With rh taking values giving 
only modest beaming factors 6>0.1 we find that the accre- 
tors in ULXs with soft components all have stellar masses 
mi < 25. I suggest also that the non-power-law ULXs obey- 
ing Laoft tx T 4 actually have black hole masses sufficiently 
high (~ IOM0) that they are SMfr-Eddington. 
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King & Puchnarewicz (2002) developed a general formalism 
for treating blackbody emission from the vicinity of a black 
hole. They allowed for geometrical collimation of this emis- 
sion, but assumed that this did not change the blackbody 
luminosity L or temperature T. This is true for example of 
radiation subject to scattering by nonrelativistic electrons. 
King & Puchnarewicz's main result (their eqn (5)) is 



L sph = 2.3 x 10 44 7T 4 



O.lkeV „lo er gS 



Here L ap h is the blackbody luminosity an observer would 
infer from the observed flux by assuming that it is isotropic, 
even though in reality it is collimated by a factor b. The 
quantity / is the ratio of the intrinsic luminosity to the 
Eddington limit Le (which can exceed unity by the fac- 
tor (1 + mm) mentioned above), p is a factor allowing 
for deviations from spherical symmetry in the source (e.g. 
that it is actually plane and inclined to the line of sight) 
and r = R/R s is the blackbody radius in units of the 
Schwarschild radius. 

The derivation of the relation ((2)1 is simple. We express 
the intrinsic (pre-collimated) blackbody luminosity as 

2 

L(xR 2 T i pocM 2 T 4 r 2 p(xL 2 T' 1 ^-, (3) 

where one writes the radius as R — rR s cx rM at the first 
step, and the mass M as M cx Le cx Ll~ x at the second. 
Solving this equation for L we find L oc T~ 4 . An observer 
assuming that the flux is isotropic with the observed value, 
rather than collimated, now infers a total blackbody lumi- 
nosity L sp h = 6 _1 L, i.e. 



I 

L sph cx T~ 4 ——, 
pbr 1 



(4) 



which gives ([2} when the proportionality constants are in- 
cluded. 

King & Puchnarewicz (2002) used (J2J) to argue that any 
source exceeding the normalization on the rhs must either be 
super-Eddington for its mass (I > 1), or emit from a region 
much smaller than the Schwarzschild radius (r < 1), or emit 
anisotropically (pb < 1). Ultrasoft quasars and some ULXs 
lie close to this regime on the L — T plane. Here, setting 
Lsph = Lsoft, we see that the L so ft cx T~ 4 correlation (jXJ) for 
ULX soft excesses is reproduced provided that 



J—=3xW- 4 . 
pur 2. 



(5) 



Observation thus strongly suggests that b oc r .If 
the power of T in fTJ were not precisely 4, e.g. the value 
n = -3.1 ± 0.5 found by Feng & Kaaret (2007), this would 
introduce a T-dependence into the relation between b and 



i.e. 6 oc T 



Since the fitted value of T 



pbr 2 



(2) 



varies by a factor < 3, while (as we shall see) the inferred 
beaming factor b varies much more, we would make only a 
small error in adopting the approximate dependence b ~ r~ 2 
here too. 

This scaling of b thus seems to be required by obser- 
vation. Theoretically, a simple argument suggests that a 
b cx r~ 2 dependence follows from the picture of ULXs pro- 
posed by Begelman et al (2006) and Poutanen et al., (2007), 
in which a wind from the accretion disc surface keeps the 
local accretion rate at the radiation pressure limit at each 
disc radius, as originally suggested by Shakura & Sunyaev 
(1973). The outflowing wind is densest near -R ap h, and has 
large optical depth both outwards along the disc plane, and 
in the vertical direction. Thus most of the disc radiation 
emitted within i? sp h must diffuse inwards by scattering, un- 
til it escapes through the central funnels parallel to the disc 
axis. The collimation results from the fact that the funnel is 
tall and thin, and has scattering walls. 

To apply the formalism leading to ([2]) we identify the 
blackbody radius R as R ~ J?s P h = rR s , with r = 27mf/4, 
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where f ~ 1. The blackbody luminosity emitted by the disc 
within -R S ph is the intrinsic luminosity L. This diffuses in- 
wards and is collimated by the funnels. 

For the beaming factor b we consider a simple cylin- 
drical funnel around the central disc axis. If the typical 
cylindrical radius of the funnel is 7?o, and its height is 
Ho 3> Ro, the half-angle over which radiation escapes is 
~ sin -1 Ro/ Ho — Ro/Ho- Then the beaming fraction b 
is simply the total fractional area of the two funnels, i.e. 
b ~ (1 — cos#) ~ Ro/2Ho- Close to the disc plane, the struc- 
ture of the central region of the disc wind (and thus the fun- 
nel radius Ro) is independent of the value of Eddington ratio 
m > 1 at large R, since all such discs have the same central 
accretion rate behaviour M(R) ~ Me(R/R s ). We expect 
that Ro ~ \R„ with A > 1, as R s sets the lengthscales in 
this region. However the funnel height Ho is sensitive to rh, 
or equivalently i? sp h = rR s , as at points far from the disc 
plane the wind flow pattern is set by rh, which is equivalent 
to saying that the large-scale flow pattern is self-similar and 
scaled by r. In particular this requires Ho ~ /-iR sp h oc /jr, 
where /i < 1, so that finally 



A 



A^ r 

2r 2 46fi 2 rh 2 



where x stands for the dimensionless combination 

I 2 



(6) 



(7) 



We see the the observational requirement ([5} implies A/yU ~ 
58, so that the funnel height is only a few percent of i? sp h, 
i.e. (i ~ few x 10 -2 . We get finally 



b ~ -—x 



(8) 



The reasoning of this paragraph assumes that rh is large 
enough that the two scales Ro and -R sp h are very differ- 
ent. The scaling of b with rh is clearly more complex for 
smaller rh. In particular, unless Ho > Ro, which requires 
m > 8.5a; -1 '' 2 , one would formally get b > 1. 



3 THE L ~ T 4 CORRELATIONS FOR 
NON-POWER-LAW ULXS 

The last Section dealt with those ULXs (the majority) for 
which soft components are seen, and obey the L S oft °t 
T -4 relation. We noted above that KP show that the re- 
maining (non-power-law) ULXs follow opposed correlations 
L ~ T 4 for the medium-energy X-rays. Here the nor- 
malization differs for each individual system. The lumi- 
nosities of the two groups differ sharply: the power-law- 
soft-excess systems have inferred luminosities permanently 
above 3 x 10 39 ergs -1 , while the non-power-law systems are 
permanently below 10 39 ergs -1 . It seems likely that these 
fainter systems correspond to sitfr-Eddington accretion on to 
black holes with masses > IOM0 . It is clear that such sys- 
tems must exist, and that there is no reason to expect the 
collimation leading to the opposite L so f t oc T -4 behaviour of 
the bright ULXs. The normalizations of the L ~ T 4 correla- 
tions are then fixed by the system inclinations and the inner 
disc radii. The latter are indeed of order a few Schwarzschild 
radii for black holes of ~ IOM0 . 




Figure 1. Beaming factor b (dashed curve) and accretor mass 
Mi (solid curve, in Mq) as functions of the Eddington ratio rh 
for ULXs. Here L40 is the inferred isotropic bolometric luminosity 
in units of 10 40 erg s — 1 and i~la dimensionless quantity given 
by equation @ 



4 MASSES, BEAMING AND EDDINGTON 
RATIOS 

We can now check whether the inferred behaviour of the 
beaming factor leads to sensible parameters for observed 
bright ULXs. Although the relation ((Sj for b was derived 
using the inferred blackbody disc emission, its geometrical 
nature and the fact that electron scattering is independent of 
photon energy makes it probable that it holds for all forms of 
ULX luminosity, and indeed even in ULXs where no black- 
body disc component has been identified, provided only that 
these correspond to super-Eddington accretion. In particu- 
lar we can use the b oc rh~ 2 scaling in considering medium- 
energy X-rays, which are generally assumed to carry most 
of the bolometric luminosity of a ULX. 

The effect of beaming is to cause an observer to infer a 
spherical luminosity 



^(l + lnm) 



(9) 



(cf Shakura & Sunyaev, 1973; Begelman et al., 2006). Elim- 
inating b using (|6} (or (J5]) gives 

36 



-1 -1 

erg s 



(10) 



(11) 



L sp h = 2.2 x lO^mim (1 + \nrh)x 
We can re-express this as 
mi _ 4500 
L40 m 2 (l + mm) 

where L40 = L sp h/10 40 ergs -1 . 

Figure 1 shows M1/L40 and b as functions of m (using 
equations [11] and [8} . We see that Eddington ratios in the 
range 8.5 < rh < 20 imply stellar masses 1M© ^ 20 for 
the accretors if the disc luminosity is < 10 41 erg s -1 , and 
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beaming factors in the range 1 > b> 0.2. Hence stellar-mass 
binaries with moderate Eddington ratios and consequently 
modest beaming provide very good candidates for explaining 
ULXs. 

We note from (|10|) that the inferred luminosity L sp h 
varies essentially only because of the sensitivity of the beam- 
ing factor b to rh. Thus the bolometric luminosity varies only 
logarithmically above Le (assuming that the Eddington ra- 
tio always exceeds unity), but is spread over a smaller or 
greater solid angle as rh increases or decreases, significantly 
altering the inferred luminosity. 

Since M = ttiMe oc {M\/rf)rh, where r], the radiative 
efficiency, is similar (~ 0.1) for black holes and neutron stars, 
I note that that some ULXs could contain neutron stars, and 
could even have lower absolute accretion rates for the same 
inferred luminosity. From Fig. 1 we see that a lOM© black 
hole with m — 15 and a neutron star with mass < 2M0 
and rh = 30 produce similar inferred luminosities, with the 
neutron-star system having an absolute accretion rate M 
lower by a factor ~ 2.5 than the black hole. The origin of this 
apparently paradox is that the latter system has a smaller 
beaming factor. (Put another way, on Fig. 1 the curves of 
constant M are hyperbolae which cross the hyperbola de- 
scribing b.) For ultrasoft ULXs with no detectable medium- 
energy X-ray component, even white dwarf accretors are 
possible, particularly since for them 77 can be enhanced over 
the pure accretion yield by nuclear burning of the accreting 
matter (cf Fabbiano et al., 2003). 



5 ULX POPULATIONS 

Population studies of ULXs have until now faced the diffi- 
culty that the beaming factor b was not determined, intro- 
ducing a spurious degree of freedom. Given the connection 
(|8|, we can now remove this. We consider a population of 
ULXs with host galaxy space density n g Mpc -3 and assume 
that each host contains N ULXs, with radiation beams ori- 
ented randomly. To be in the beam of one such object one 
has to search through ~ 1/Nb galaxies, i.e. a space volume 
~ l/n g Nb. The nearest observed ULX is thus at a distance 



D„ 



AnrigNb 



1/3 



0.7(n 9 iV) 



-1/3 • 2/3 



Mpc, 



(12) 



where rhi = m/10. The apparent luminosity of the ULX is 



L sph = 2.2 x 10 39 m*rh 2 ergs" 1 



(13) 



where m„ = Mi/IOM©, giving a maximum apparent bolo- 
metric flux 



Fmax — 



^ = 4.0xl0- 11 m»m? /3 (n 9 iV) 



2/3 -1 -2 
erg s cm 



(14) 



These relations, together with the results of the previous 
Section, agree with the fact that ULXs of apparent lumi- 
nosity few x 10 39 - 10 41 erg s 1 are observed in the Local 
Group, and suggest that the typical intrinsic number N per 
host galaxy is at most a few. This is in line with estimates 
of the numbers of high-mass X-ray binaries in phases of 
rapid mass transfer on thermal or nuclear timescales (King 
et al., 2001; Rappaport et al., 2005), suggesting that these 
systems offer good candidates for explaining most if not all 



ULXs. Ultimately one needs a population synthesis calcula- 
tion to verify that this picture produces the right numbers 
of systems with the required moderate Eddington ratios to 
produce the nearby ULXs. 



6 PSEUDOBLAZARS? 

It is unclear to what value of rh one may safely extrapolate 
the b oc rh~ 2 dependence inferred here. This is an interest- 
ing question, as we know (cf Begelman et al, 2006; King 
& Begelman, 1999) that the well-studied object SS433 has 
rh ~ 3000 — 10 4 . Such values are typical for both thermal- 
timescale and nuclear-timescale mass transfer from massive 
donor stars (Rappaport et al., 2005). 

From the work of the previous Section, now scaling rh as 
rh = 10 4 rfi4, the nearest such object would be at a distance 



AnrigNb 



1/3 



6607V" 1/3 m4 /3 Mpc, 



(15) 



where I have taken n g ~ 0.02 Mpc -3 as appropriate for L* 
galaxies. The apparent isotropic luminosity of such an object 
would be 



L sp h = 2.2 x 10 4j m*rh 2 ergs 1 . 



(16) 



Hence in distance and apparent luminosity the object would 
appear as an AGN. However, unlike a genuine AGN, there 
is no requirement that it should lie precisely in the nucleus 
of the host galaxy. 

A possible candidate for such an object is the BL Lac 
system PKS 1413+135 (Perlman et al., 2002). With redshift 
z = 0.24671 it has distance D ~ 1000 Mpc and isotropic 
luminosity ~ 10 44 ergs -1 , but lies at 13 ± 4 mas from the 
centre of the host galaxy. 



7 DISCUSSION 

The work of this paper suggests that the beaming factor in 
super-Eddington accretion varies as b oc m~ 2 . This seems 
to be required by observations of the L so ft — T correlation, 
and is reasonable on general geometrical grounds. The ex- 
istence of this scaling means that observable properties of 
ULXs depend essentially only on the Eddington ratio rh. If 
this conclusion is valid, this removes the spurious degree of 
freedom allowing one to choose b independently of rh which 
has made systematic parameter estimates difficult in the 
past (e.g. King, 2008, where these two quantities are not 
connected). 

It appears that most ULXs correspond to stellar mass 
systems accreting at Eddington ratios of order 10 - 30, with 
corresponding beaming factors 6>0.1. High-mass X-ray bi- 
naries containing black holes or neutron stars are excel- 
lent candidates, although population synthesis studies are 
needed to check this. The scaling inferred here suggests that 
ULXs vary mainly because the beaming factor is sensitive 
to the Eddington ratio. If the scaling can be extrapolated 
to the Eddington ratios found in SS433, some objects cur- 
rently identified as AGN at modest redshifts might actually 
be ULXs. This may explain cases where the AGN does not 
coincide with the centre of the host galaxy. 
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